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ABSTRACT 

We calculate the column density distribution of damped Lya systems, modeled 
as spherical isothermal gaseous halos ionized by the external cosmic background. The 
effects of self-shielding introduce a hump in this distribution, at a column density Nhi ~ 
1.6 X W^'^X-^ cm , where X is the neutral fraction at the radius where self-shielding 
starts being important. The most recent compilation of the column density distribution, 
by Storrie-Lombardi & Wolfe, shows marginal evidence for the detection of this feature 
due to self-shielding, suggesting a value X ~ 10~^. Assuming a photoionization rate 
r ~ 10~^^s~^ from the external ionizing background, the radius where self-shielding 
occurs is inferred to be ~ 3.8 kpc. If damped Lya systems consist of a clumpy medium, 
this should be interpreted as the typical size of the gas clumps in the region where 
they become self-shielding. Clumps of this size with typical column densities Nh ~ 
3 X 10^'^ cm~^ would be gravitationally confined at the characteristic photoionization 
temperature ~ 10^ K if they do not contain dark matter. Since this size is similar to 
the overall radius of damped Lya systems in Cold Dark Matter models, where all halos 
are assumed to contain similar gas clouds producing damped absorbers, this suggests 
that the gas in damped absorbers is in fact not highly clumped. 

Subject headings: galaxy: formation - quasars: absorption lines 



1. INTRODUCTION 

Damped Lya systems (hereafter, DLAs) contain most of the atomic hydrogen in the universe 
(e.g., Wolfe et al. 1995; Storrie-Lombardi & Wolfe 2000). They are therefore fundamental objects 
to understand galaxy formation: in the process of forming galaxies, the matter that is originally 
in the low-density, ionized intergalactic medium must invariably become first atomic, before it can 
form molecular clouds and stars. 
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The nature and geometrical shape of the DLAs is still a subject of debate. The two alternative 
hypotheses that are being discussed at present and confronted with the observations are as follows: 
(a) They arc rotating disks, which may or may not be forming stars; (b) They are approximately 
spherical halos, where global rotation contributes little to the support against gravity. The truth 
may be a combination of these two possibilities: DLAs might be described as flattened halos or 
thick disks, with a range of the ratio of rotation velocity to velocity dispersion. 

A fundamental property of DLAs is that their associated metal lines reveal the presence of 
multiple absorbers in most of the systems, which are thought to arise from clumps of gas. The mul- 
tiple absorption lines appear over velocity intervals of 30 to 300kms~^. The statistical properties 
of the multiple absorption lines have been used to attempt to distinguish between the disk and halo 
models, but both alternatives seem to be consistent with observations (Prochaska & Wolfe 1997, 
1998; Haehnelt, Steinmetz, & Ranch 1998, 2000; McDonald & Miralda-Escude 1999). Although 
the observations reveal that there are usually one or a few clumps intersected along a random line 
of sight, with characteristic column densities of ~ 10^°cm~^, the size of the clumps is not well 
known. Constraints in a few systems from double lines of sight in lensed QSO's, and from model 
calculations of the ionization parameter, indicate sizes or structure in the range of 20 pc to 1 kpc 
(Rauch, Sargent, & Barlow 1999; Lopez et al. 1999). Models of structure formation predict that 
the overall size of DLAs is of the order of 1 to 10 kpc (e.g., Katz et al. 1996; Gardner et al. 1997; 
McDonald &; Miralda-Escude 1999). If the clump size is not much smaller than the size of DLAs, 
then the absorption lines would actually be arising from mild fluctuations in density and velocity 
in a turbulent medium; whereas if the clump sizes are much smaller, they would be real separate 
entities with a large overdensity relative to an interclump medium. 

This paper discusses the effects of self-shielding on the form of the column density distribution 
of damped Lya systems, and the dependence of this distribution on the size of the gas clumps. 
As the gas becomes self-shielding against the external cosmic ionizing background, there is a rapid 
transition from ionized to fully atomic gas, and this causes a feature in the column density distri- 
bution. If X is the neutral fraction of the gas seen at a column density Nhi = 1.6 x lO^^cm"^ 
(where the optical depth to photons at the hydrogen ionization edge is equal to one), a hump in the 
column density distribution will occur at Nhi ~ 1-6 x lO^^X'^ cm . Therefore, a measurement 
of X from this feature in the column density distribution, combined with an independent estimate 
of the ionizing background intensity, can be used to infer the size of the gas clumps. 

Observationally, compared with the expected number extropolated from low column density 
distributions, an excess of absorption systems with Nhi ^ 2 x lO^*' cm~^ is found by Lanzetta et 
al. (1991). The recent compilation of the column density distribution of known DLAs by Storrie- 
Lombardi & Wolfe (2000) shows more clearly a hump around Nhi ~ 2 x 10^*^ cm~^. Theoretically, 
based on an approximate calculation of the self-shielding effect, Murakami &; Ikeuchi (1990) showed 
that a flat part appears in the column density distribution. With a simplified consideration on the 
ionizing flux transfer, Petitjean, Bergeron, & Puget (1992) investigated models of self-gravitating, 
photoionized, spherical gaseous cloud and explained the flattening of the column density distribution 
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at Nhi ~ 2 X 10^° cm~^ as the effect of self-shielding. Similar calculations of the expected column 
density distribution due to self-shielding have been done by Corbelli, Salpeter, &; Bandiera (2001) 
for plane-parallel geometry. 

Here, we will consider a spherical geometry with a singular isothermal profile to perform a 
self-shielding calculation, and focus on the application of inferring the size of the gas clumps in 
damped Lyct systems. 

2. METHOD 

2.1. Model for Gaseous Halos 

We model a DLA as a spherical cloud of gas with a singular isothermal profile, pg = A/r^, 
with a cutoff at the virial radius, where pg is the gas density. The model depends only on the 
constant A, which can be recast in terms of the mass of an associated virializcd dark matter 
halo. We will be presenting results for dark matter halo masses of 10^, lO^'', 10^^, and 10^^ M©, 
assuming that 5% of the mass of the halo is in the form of gas. For the cosmological model 
with CI = 1, Hq = 70kms~^ Mpc~^, and at z = 3, the virial radii corresponding to these halo 
masses are Tvir = 5.2,11.1,23.9, and 51.5 kpc, and the halo circular velocities are Vc = 29,62,134, 
and 289kms~^, respectively (e.g., Padmanabhan 1993). In Cold Dark Matter models of structure 
formation, the damped Lya systems should be located in halos over this range of masses (e.g., 
Gardner et al. 1997). The constant A is simply determined by the condition that 5% of the halo 
mass is equal to the total gas mass within radius ryir- The model ignores dumpiness of the gas, 
but as discussed later the results depend only on the assumed radius at a given density, so the halo 
radius will essentially play the role of the radius of a clump. 

2.2. The Background Spectrum And The Ionization Profiles 

The other quantity that determines the profile of neutral hydrogen is the intensity of the 
external ionizing background. We assume an intensity of J,^ = 3 x 10^^^ ergcm^^ s^^ Hz^"*^ sr~^, 
constant at all frequencies between the Hi Lyman limit, ul, and the Hell Lyman limit, 41^^. 
The intensity is set to zero at frequencies above Av^. This is a sufficiently good approximation 
for our purpose to the background spectrum expected from QSOs, after absorption by the Lya 
forest is taken into account (e.g., Haardt &; Madau 1996). The reason is that Hen becomes self- 
shielding at a larger radius in the gaseous halo than Hi (owing to the faster Hell recombination 
rate and the lower intensity of the cosmic background at the Hell Lyman limit, compared to 
hydrogen; e.g., Miralda-Escude & Ostriker 1990), and therefore the photons at frequency above Aul 
in the cosmic background are absorbed before reaching the region where Hi opacity is important. 
With this spectrum, our adopted value of the intensity at the Hi Lyman limit corresponds to a 
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Hi photoionization rate T = 1.25 x 10 s . 

In our calculation, we will set the total hydrogen number density to uh = 0.76pg/mH (where 
niH is the mass of the hydrogen atom), and the electron number density to Ug = 0.82(1 — 
XHi)pg/n^Hj where xhi = nni/riH is the hydrogen neutral fraction. The factors 0.76 and 0.82 
assume a helium abundance of y = 0.24 by weight, and that all the helium is in the form of He I 
and Hell with the same neutral fraction as hydrogen. Other than these factors, the presence of 
helium is ignored in our calculation. This is a good enough approximation owing to the small 
helium abundance by number, plus the fact that every photoionization of Hei produces a photon 
following recombination of He ii which ionizes hydrogen. 

2.3. Self-shielding Calculation 

We compute the self-shielded H i density profile using an iterative procedure similar to the one 
described by Tajiri &; Umemura (1998). Photoionization equilibrium demands that at each radius, 
we have 

xm\r) / d\l \ a^du = [1 - XHi{r)\ ^ aB{T) , (1) 

Jatt Jul ''^H 

where xhi{t) is the hydrogen neutral fraction at radius r, Ivfi is the cosmic background intensity, 
is the photoionization optical depth from radius r to infinity along a given direction, ai, is the 
photoionization cross section of hydrogen, and Q!b(T) is the case B recombination coefficient at 
temperature T (we use the on-the-spot approximation, see Osterbrock 1989). We fix the temper- 
ature to 2 X 10*^ K, for which the recombination coefficient is as = 1.41 x 10~^^cm'^s~^ (Spitzer 
1978; Verner &; Ferland 1996). One can compute self-consistently the gas temperature assuming 
thermal equilibrium under photoionization conditions, but in practice the gas is likely to be fre- 
quently shock-heated in any model of moving gas clumps in a halo. To numerically evaluate the 
integral in equation (1), we compute the integral over frequency at many values of the optical depth 
and obtain a fit to this function. The integral over solid angle is reduced to one dimension in the 
spherical case, and is performed at each radius. We initiate XHi{r) by assuming the system to be 
optically thin. At each iteration, is recomputed for each radius and polar angle, and new values 
of Xhi are found. This iteration is repeated until xhi has converged with a fractional error less 
than 10-^. 

3. RESULTS 

The neutral fraction profile as a function of radius, in kpc, is shown in the left panel of Figure 1 
for the four halo masses. As expected, the ionizing photons penetrate deeper in halos of lower mass. 
For a power-law density profile, self-similarity implies that the shape of the xhi profile depends on 
the gas density and the background intensity only through one parameter, X, which is the neutral 
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Fig. 1. — The neutral fraction profile as a function of radius for four halo masses. As explained 
in the text, these xhi profiles depend only on the variable X, which is the neutral fraction at the 
self-shielding radius rss- The radius is in physical unit for the left panel and is normalized by Vss 
for the right panel. 

fraction at the radius Vgg where the cloud becomes sclf-shiclding: X = XHi,oiTss)- Here, rgs is 
defined as the impact parameter where the column density is Nffiflirss) = — ^-^^ ^ ^^^^ cm~^ 
with both xhi,o and Nnifl computed without including the self-shielding correction. Using the 
approximation X <C 1 and ignoring also any external cutoff in the isothermal gas density profile, 
one finds 

y 0.767rrrss 

Thus, the shape of the neutral fraction profile depends only on the product Fr^^. It is easily shown 
that the self-shielding radius Vss is proportional to (A^/F)^/^ oc M^/^/F-*^/^. For the four halo 
masses we use, (lO^, 10^°, 10", and lO^^M©), and F = 1.25 x IQ-^^g-i^ ^ 1.3,3.6,9.9, 

and 27.6 kpc, and X = 1.8 x 10"■^ 1.1 x 10-3,6.3 x 10"^, and 3.8 x lO""^, respectively. The right 
panel of Figure 1 shows the neutral fraction profiles for these four cases as a function of r/rgs- 

This general shape of the neutral fraction profiles of self-shielding clouds ionized by the cosmic 
background can be simply understood as follows: the cloud becomes self-shielding at the radius 
where the total rate of recombinations balances the rate at which photons come into the cloud. 
Therefore, once a radius ~ r^^ is reached where most photons with frequency near z/j;, have been 
absorbed, the cloud becomes mostly neutral over a small range in radius, because the rate of 
recombinations required to balance the absorption of all the higher frequency photons is not much 
higher. The effect has been demonstrated before by Murakami & Ikeuchi (1990) and Petitjean, 
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Bergeron, & Puget (1992), and is analogous to the sharp edges of Hi disks produced by self- 
shielding (Corbelli &: Salpeter 1993). Since 04,^^ ~ a^^/QA, we conclude that for X > 1/64, some 
of the highest frequency photons will penetrate into the atomic region of the cloud (causing a 
residual ionization inside the atomic region), while for X < 1/64, all the photons have mean free 
paths small compared to rgg once the cloud is mostly atomic, so the intensity of ionizing photons 
suddenly drops to zero and there is a sharp transition to a completely atomic gas. 

Figure 2 shows the predicted distribution of Hi column densities, /{Nhi), when random lines 
of sight intersect the four model gas clouds (obtained by requiring /{Nhj) dNni oc rdr, where 
Nhi{t) is found by integrating the hydrogen density times the neutral fraction in Fig. 1 over lines 
of sight at impact parameter r). The transition to atomic gas causes a feature in the distribution, 
with a hump appearing at column density Nhi ~ Nnifl . 




10" 1018 1019 IQ20 IQZl IQZZ 

Nh, (cm-2) 

Fig. 2. — The probability distribution of the neutral hydrogen column density for different halo 
masses. Each curve is normalized so that the total probability of intersecting the halo within the 
virial radius is unity. Data points are from Storrie-Lombardi & Wolfe (2000) with a vertical shift. 

Superposed on our curves, we reproduce a recent compilation of the column density distribution 
of known damped Lya systems, from Storrie-Lombardi & Wolfe (2000). The points seem to fit best 
the curve for the 10^^ M© halo (notice that the points can be arbitrarily shifted vertically). In 
particular, the point at the lowest column density suggests a shallower slope of the distribution. 
This is a marginal indication that the feature predicted by the self-shielding effect may have been 
detected. If so, the value of X that is inferred is X ~ 1.1 x 10~^, which for our adopted value of the 
photoionization rate F, corresponds to r^s ~ 3.6 kpc. We caution, however, that the predicted shape 
depends on the assumed total gas density profile, which we have taken to be singular isothermal. 
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4. DISCUSSION 

We have shown that self-shielding should cause a sharp transition from an ionized to an atomic 
medium in any gas cloud photoionized by the external cosmic background. The expected absence 
of photons above frequency Au^ due to He ii absorption makes that transition even sharper once the 
neutral fraction has been increased by a factor 64 by sclf-shiclding. This produces a hump in the 
column density distribution. The column density at which the hump occurs measures the quantity 

rrss. 

Most damped Lya systems are likely to be photoionized by the external cosmic background. In 
fact, since their rate of occurrence is about one fourth per unit redshift at z ~ 3 (Storrie-Lombardi, 
McMahon, &; Irwin 1996), and the mean free path of ionizing photons is Az ~ 1 (Miralda-Escude 
Sz Ostriker 1990), it is easily shown that if all the sources of the cosmic ionizing background were 
embedded inside damped Lya systems, then the local source in any such system would contribute 
about the same flux as the external background on average. It is much more likely that any 
sources associated with typical damped Lya systems do not contribute significantly to the cosmic 
background. 

The results for the column density distribution presented by Storrie-Lombardi & Wolfe (2000) 
suggest a possible detection of the self-shielding effect, which would imply a radius rgs — 3.6 kpc 
for r = 1.25 X 10~^^s^^. This represents a measurement of the size of the clumps in damped 
Lya systems. Even though our calculation is for a spherical, isolated cloud, a system of randomly 
located clumps would produce a similar column density distribution depending on the radius of the 
clumps, as long as each clump is illuminated by the ionizing background with an optical depth < 1 
along a large fraction of directions. This should be correct since the typical number of intersected 
clumps identified in the metal absorption lines is not very large. 

The typical size of damped Lya systems inferred from CDM models where all halos with 
Vc > 40kms~^ give rise to the absorption systems is 1 to 10 kpc (e.g., Katz et al. 1996; Gardner et 
al. 1997; McDonald & Miralda-Escude 1999). Hydrodynamic cosmological simulations show that, 
at redshift 2; ~ 2 — 4, for halos more massive than ~ IO^^Mq, the projected distance of DLAs to 
the center of the nearest galaxy is around 10 kpc (Gardner et al. 2001). Since the clump size we 
infer is about the same, this suggests that the multiple metal absorption lines do not arise from 
highly overdense clumps, but from mild density fluctuations in gaseous halos. If there were highly 
overdense clumps, their small self-shielding radii would produce a column density distribution with 
a less pronounced hump at lower column density than the curves in Figure 2. The appearance of 
individual absorption lines in the spectra of Mg 11 and other metal lines might be deceiving, and 
could be due to a continuous medium, where the Mg 11 density is highly sensitive to the gas density 
due to photoionization and self-shielding effects. 

Nevertheless, a number of caveats must be borne in mind in interpreting the observed column 
density distribution as a measurement of Frgs. In reality, there should be a distribution of clump 
sizes, and differences from a simple spherical geometry would introduce additional variations, which 
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would smooth the shape of the distribution predicted in Figure 2. Furthermore, the density profile 

of the gas cloud can be different from what we assume here. A fully three-dimensional calculation 
of self-shielding in a hydrodynamic simulation of galaxy formation would be highly desirable for a 
more robust interpretation of the observations of damped Lya systems. 

Finally, we mention that once the size of the clumps is known, their column densities can be 
used to infer their gas mass, and the temperature required to have a cloud in hydrostatic equilibrium 
(see Corbelli, Salpeter, & Bandiera 2001). Using a typical column density A^^f ~ 3 x 10^*^ cm~^ 
for a clump, and a radius Vgs — 3.6 kpc, a gas mass of Mg ~ 1.6 x 10^ is inferred for the 
self-shielded region. In the absence of dark matter, the temperature in hydrostatic equilibrium is 
T ~ GMymH/krss ^ 2.3 x lO'^K. Since this is near the expected temperature in photoionization 
equilibrium, this would imply that the gas clumps in damped systems do not contain a lot of dark 
matter. However, this conclusion is altered if the gas temperature is higher due to shock heating 
or turbulent motions are important. 

This work was supported in part by NSF grant NSF-0098515. 
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